Majorana edge modes of topological exciton condensate with superconductors 
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I study the edge states of the topological exciton condensate formed by Coulomb interaction 
between two parallel surfaces of a strong topological insulator. When the condensate is contacted 
by superconductors with a n phase shift across the two surfaces, a pair of counter-propagating 
Majorana modes close the gap at the boundary. I propose a nano-structured system of topological 
insulators and superconductors to realize unpaired Majorana fermions. The Majorana signal can 
be used to detect the formation of the topological exciton condensate. The relevant experimental 
signatures as well as implications for related systems are discussed. 



Introduction. — Topological states of matter exhibit a 
novel form of quantum order, leading to precise quanti- 
zation of certain physical quantities [T]. Important ex- 
amples of such topologically ordered states are provided 
by the hierarchy of quantum hall states [2H1] and time- 
reversal invariant topological insulators (TIs) [IHZ]- The 
combination of topological order with the conventional 
order arising from broken symmetries gives rise to new 
degrees of freedom, such as protected chiral or helical 
surface states and Majorana modes, that are otherwise 
not realized in nature. The study of these states is cru- 
cial to our understanding of the collective properties of 
matter and could be important for applications ranging 
from high-precision metrology to quantum computation. 

Recently, a topological exciton condensate (TEC) was 
predicted by a group including the author to exist when 
a thin film of strong TI is gated so that electron and 
hole gases separately form on opposite surfaces due to 
Coulomb interaction |8l [9]. In the exciton condensate 
state, the system acquires a coherent tunneling ampli- 
tude between the two surfaces even though there is no 
direct electronic tunneling between them. It was shown 
that the topological nature of the TEC gives rise to frac- 
tionally charged vortices of the condensate. A similar 
system based on graphene has also beed proposed and 
studied [TUHT2]; however, due to the higher multiplic- 
ity of Dirac cones, the exciton condensate in graphene is 
topologically trivial. 

In this Letter, I study the structure of edge states 
when the TEC is interfaced with other ordered states 
and propose a novel way to detect the TEC experimen- 
tally. Remarkably, I show that under certain conditions 
a (nonchiral) Majorana channel exists at the boundary 
of the TEC with a superconductor that spans the two 
surfaces. In the lab, it might be better to substitute 
the thin film with a dielectric wafer sandwiched between 
two TIs. This design offers greater control on the dielec- 
tric screening of the insulating spacer [T3] . Fig. fl] shows 
the schematics of the proposed setup for both a thin-film 
and a sandwich structure. As the superconducting phase 
shift across the Josephson junction is cycled through 27r, 
the Majorana channel opens once, resulting in two un- 
paired Majorana modes moving along the channel in op- 



FIG. 1. Majorana fermions. The TEC is attached to su- 
perconducting leads and contacts with a phase difference Sip 
between the top and bottom surfaces, (a) top view, and the 
proposed setup in (b) thin film TI and (c) sandwich struc- 
tures, where a dielectric wafer is sandwiched between two TI 
slabs. By tuning Sip in the middle contact, a pair of counter- 
propagating Majorana modes (arrows) form at the boundary 
with the TEC, resulting in two unpaired Majorana fermions 
(dots) at the intersection. 



posite directions. These Majorana modes exhibit non- 
Abelian fractional statistics [2] and may be detected us- 
ing a suite of experimental techniques that are currently 
used or being developed p^Hl8] . This prediction can be 
exploited either to realize and manipulate unpaired Ma- 
jorana fermions or, more significantly, to detect the TEC 
itself. 

The existence of the Majorana channel may be un- 
derstood in terms of the Majorana edge modes on an iso- 
lated surface of a TI, which are found in the following two 
cases. First, there is a pair of counter-propagating Majo- 
rana edge modes at a domain wall separating (proximity- 
induced) superconducting regions with a phase difference 
TT [19] . As shown in Fig. l2](a) by folding this domain wall 
across the side surface to the opposite surface these Ma- 
jorana edge modes map to the ones described above. Sec- 
ond, there is a chiral Majorana edge state at the bound- 
ary of a superconductor on a single surface and a mag- 
netic region dominated by the Zeeman energy. As shown 
in Fig. |2jb), by folding the magnetic region as before, 
we obtain a pair of counter-propagating Majorana edge 
modes in each surface living at the boundary between su- 




FIG. 2. Constructing the Majorana edge modes by bending 
a single surface, starting from (left) a domain wall in the 
single-layer superconductor with (^1 — 792 ~ tt, and (right) the 
boundary between superconducting and magnetic regions. 



perconducting and antiparallel magnetic regions on the 
two surfaces. The mass term associated with the antipar- 
allel magnetic order, indicated by M^ in Eq. (fTl) below, 
anticommutes with the TEC mass term, i.e. they add in 
squares to give the square of the energy gap in a region 
with both orders. Therefore, we may adiabatically switch 
the TEC on and the magnetic region off without closing 
the bulk gap. Since this process couples the two surfaces 
through exciton tunneling terms, the pair of Majorana 
modes will survive only if protected by a symmetry. I 
will show that this is precisely what happens when the 
phase of the superconductor between the two surfaces 
is adjusted to tt and the symmetry is closely related to 
time-reversal (T) symmetry. 

Mass terms. — The Hamiltonian describing the surface 
states of the system is H — ^j.-0^(i")ft-o'0(r) where 



ho = Tz{v(T • p + V") + \m\Tx:e 



«XT^ 



M, 



(1) 



and the spinor tp'^ = e*'^'^^/'*(V'it, V'l^i '/'2t; V'2^) has four 
components indexed by surface (a = 1,2) and spin (t,4) 
labels, with v the Fermi velocity, t and cr Pauli matrices 
acting on the surface and spin index, respectively, and 
p — —i{dj;,dy) the momentum operator. Here, \m\e^^ is 
the complex exciton order parameter, V is a symmetric 
bias, and I have included the magnetic order parameters 
M = Mpaz + MaTzCTz for the parallel (Mp) and antipar- 
allel {Ma) components along the surface normal. 

To include the superconducting order parameter, 
we pass to the Nambu spinor ^^ — {%l)'^ ^iip^ay) 
and the Bogoliubov-de Gennes Hamiltonian h = 

A -(Jyhlay 



The superconducting order parame- 
ter in each surface has the form |AQ,|e"''°, so that 



A== |A|e''^A^-e'*'^^'-/2, 



(2) 



where |A| = a/JA^A^, A = ^/\K^JK^\, the total phase 
tp = 2('^i + f2) and the phase shift Sip = Lpi — (/92, and 
77 are Pauli matrices acting on the Nambu index. I set 
A = 1, for now, and consider other values later. We may 
write the full Hamiltonian as h = U\hs + M)U, with 



which shows only the combination Sip — 2x is physically 
significant. 

These eight mass terms leave out two additional su- 
perconducting mass terms ~ VxTy'^z and rjyTya^ that 
mediate pairing between the two surfaces. They may 
be generated either by interlayer proximity effect via the 
superconductor or by an interlayer pairing interaction. I 
assume these processes are negligible. 

Symmetries. — The T symmetry is given by the antiuni- 
tary operator Q = idyK with K the complex conjugation 
and is broken for a general value bip — 2x ^ Ini:., n € Z. 
Physically, these values amount to having an inplane flux. 
For 8p ~2x = (2n -I- l)7r, anticommutes with the last 
term in Eq. (pi) . Note that r\z also commutes with all the 
terms except the last term with which it anticommutes. 
Therefore at these special values, a new symmetry is ob- 
tained given by the antiunitary operator T = 77^ 9. This 
will be important for discussing the edge modes. 

Majorana edge modes. — The edge modes for various 
phase boundaries can be understood in terms of the al- 
gebraic relations of the corresponding mass terms; if they 
anticommute, the gap will never close at the boundary 
and hence there are no gapless edge modes; if they com- 
mute, the gap closes at the boundary and a gapless edge 
channel opens. To see this, let us consider a boundary 
along the y axis with |r7T,| — |A| < for a; — 7> 00 and > 
for X — > —00 and set M = 0. For V^ = 0, wc find the edge 
modes 

$(p ) = (7g'P!/!/g/°'(l™('')|-|A(r')sin(5(^/2-x)|]dr/i)^ r^\ 

where Tj^cr^^o = s77^crx*o = -*o, s = sgn[sin((5v3/2-x)], 
and C is a normalization factor. This eigenvalue problem 
has two solutions, which can be chosen as eigenstates 
VzTz'^y'^o = i'J'o 1 since rjzTzdy commutes with both 
TyGx and rfxCx. Projecting to the subspace spanned by 
^J on which Pauli matrices p act, the Hamiltonian takes 
the form 



hs\o = vpzPy +Scos{Sp/2-x)Px, 



(5) 



where the overlap (5 = / ^g ^{x)rix'iQ{x)\A{x)\dx. The 
energy of the bound states is then found to be 



E{Py 



± 



v^pI 



- 6'^cos'^{Sp/2 — x)- 



(6) 



hs = r]zTz{v(T-p+V) + \m\rizTx + \A\r]xe 



i{5ip/2-x)rizr^ 



, (3) 



When 6p — 2x = (2n + l)7r, the edge states form 
a pair of counter-propagating gapless Majorana modes 
along the edge. The degeneracy a,t py = is protected 
by T symmetry. The two solutions ^q are T-partners, 
i.e. T^g = ivf^. Since T^ — —1, Kramers theorem 
applies and ensures that a pair of zero-energy states sur- 
vive in the presence of potentials that do not break T 
symmetry. Since T is fundamentally connected to @, 
it seems this protection should be equally robust as the 
one under T symmetry. As a direct consequence, since 
TryzTzT"^ = r]zTz, a non-zero y-term will not split the 




FIG. 3. Unpaired Majorana fermion. Starting from (a) a 
droplet of TEC, one (b) condenses a pie of the superconduct- 
ing state with the phase 5(^ winding from to the pie angle. 
At (c) Sip — 2x = 1^ and a pair of Majorana modes form along 
the radius. When (d) the droplet is filled with the supercon- 
ductor with a 27r vortex in Stp an unpaired Majorana fermion 
sits at the vortex core and one at the boundary. 



degeneracy. A A 7^ 1 introduces an additional term 
~ Tjy in Eq. (pi) for 5(p — 2x — (2n -|- 1)tt. Again, since 
ijy — TrjyT^^, the gapless edge states persist. 

Unpaired Majorana fermions. — The above discussion 
suggests that an unpaired Majorana fermion must exist 
when Sip contains a vortex of winding 27r. Imagine start- 
ing with a circular droplet of TEC, as in Fig.[3]and slowly 
opening a pie with the superconducting order parameter 
where Sip varies from up to the pie angle. As the an- 
gle ranges from to 27r, at some point along the process 
the condition 6ip — 2x = t^ is satisfied and a channel of 
counter-propagating Majorana modes open along the ra- 
dius, thereby localizing a Majorana fermion at the center 
of the droplet and another at the boundary. Finally, the 
pie closes onto itself, creating a vortex of winding 27r in 
6ip and leaving a single Majorana fermion at the vortex 
core and its partner at the boundary. 

In Hamiltonian pi), since we have mapped to tp — 0, 
the vortex in 5ip consists of ±7r windings of the super- 
conducting phase for each surface. Therefore, this sym- 
metric description introduces branch cuts in the order 
parameter. A nonsingular description is obtained if (p 
also contains a ±7r winding, so that there is a full quan- 
tum vortex in one surface and none in the other. In other 
words, the vortex in Scf) appears as a magnetic monopole 
with flux h/2e penetrating one surface and then spread- 
ing between the planes. The single Majorana mode then 
lives in the surface containing the vortex [19]. Our dis- 
cussion shows that this Majorana bound state is stable 
in general when there is a vortex in 6ip> regardless of how 
it is distributed between the two surfaces. 

We can check this claim explicitly, e.g. for the junction 
geometry in Fig. Ha) where the exciton order parameter 
is \m\Txe^^'^' for a: > 0, |y| < V^\x\ and zero otherwise, 
and the superconducting order parameter is \A\e'^^'^'^^^'^ 
for x < 0, \y\ > v3|a::| and zero otherwise with 6ip = 6(fa 
for y > 0, and Sipt, for y < 0. Assuming the junctions 
have vanishing width and setting y = 0, we see from 
Eq. (p]) that a bound state exists when there is a kink in 
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FIG. 4. (a) Unpaired Majorana fermion (dot) at the junction 
of two superconducting contacts and the exciton condensate 
(top view), (b) The topological phase diagram where regions 
with orthogonal unpaired Majorana fermions are shaded and 
hashed. The solid and open dots mark the values at which 
an analytical solution for the Majorana bound state is given 
in Eq. Q. 



the px mass term, i.e. when cos{5ipa/2 — x) cos{6ipi)/2 ~ 
X) < 0. For Sipa — 2x = tt/S and Sipi, ~2x — 57r/3, shown 
by the dot in Fig. WCb), and |A| — \m\ an analytical 
solution for the Majorana bound state, hgT = 0, may be 
found as 



p — (^g-|™|n-r/t> 



7, 



(7) 



where 7 = 4=(^(j"— ic^g ) and n is a unit vector bisecting 
each region as shown in Fig. Ka) . There is only one solu- 
tion and the sign c = sgn[cos((5(/9a/2 — x)] — +1 is chosen 
by the orientation of the superconducting regions and the 
exciton condensate. Finally, non-zero V will not remove 
this zero-energy state since the spectrum is particle-hole 
symmetric, flhsfl~^ = —hg with Q. = rjyUyK . 

The bound state is stable as 5ipa and 5ipb change so 
long as the gap along the boundaries does not close. The 
Hamiltonian ([3]) is 47r-periodic in Sip — 2x- This is a man- 
ifestation of the fractional Josephson effect due to Majo- 
rana modes [501 121] ■ Note that when \5'pia — 5pb\ = 27r 
there will be a domain wall along the boundary of the 
two superconductors. This leads to two nonchiral Ma- 
jorana channels along that edge. While this does not 
change the parity of the unpaired Majorana fermions, it 
does change the superposition of ^q in 7 to an orthogo- 
nal combination. For example, c changes from +1 to —1 
when 5ipa — 2x — 77r/3 and 5ipa — 2% = IItt/S, shown by 
an open dot in Fig. l4l^b). Combining these facts, we find 
the topological phase diagram shown in Fig.Qb). 

Discussion. — Although expected from our general ar- 
guments, it is still remarkable that we find a single Majo- 
rana bound state for an 8 x 8 underlying Hamiltonian, h^. 
This is the same dimensionality as the low-energy Hamil- 
tonian of superconducting state of spinless fermions on 
a honeycomb lattice (such as graphene), where the role 
of spin and exciton order parameter is played by the 
sublattice pseudospin and the Kekule distortion, respec- 



tively. Previous studies found an even number of Majo- 
rana fermions in such a system [22l [23] . In our case, the 
superconducting order parameter with phase difference 
dip across the two surfaces translates to a superconduct- 
ing state with intra- valley pairing and phase difference S^p 
across the two valleys. The pairs in this superconducting 
state of graphene have finite momentum. A related sys- 
tem is the new family of iron pnictide high-temperature 
superconductors where there is evidence for s* pairing. 
Indeed, a recent study [24! reached similar conclusions in 
a lattice model with an 8 x 8 Hamiltonian. However, it 
seems difficult practically to engineer a vortex configura- 
tion in these systems, since the phase shifts that need to 
be tuned occur in momentum space. 

These findings can be utilized in two ways. First, they 
can be used to create and manipulate Majorana fermions. 
The protocols for doing so are similar to those proposed 
by Fu and Kane P^. However, in regard of the increased 
complexity of the device designs, this is not a particu- 
larly advantageous application. Second, and more signif- 
icantly, the detection of the Majorana modes can be used 
as a proxy for the TEC. This opens the way to new tech- 
niques p!5Hl8] to signal the formation of the TEC, which 
are not possible for other schemes of creating an exciton 
condensate, such as quantum Hall bilayers, graphene, or 
semiconductor quantum wells. I will briefly outline the 
procedure, discussing separately the case of a thin film 
and the sandwich structures of Fig. [l] In all these cases, 
controlling the phase shift Sip — 2x is vital. This can 
be achieved by, say, an in-plane field, along the intended 
Majorana channel. 

In a thin film, when the two opposite surfaces are so 
close that there is direct electronic tunneling between 
them, the bulk TI is lost. There will be no Majorana sig- 
nal in this case as its existence relies on there being helical 
surface states of the bulk TI. (Note that in the thin film 
the superconducting contacts must be placed outside the 
tunneling region.) The TEC is expected to form when 
the two surfaces are far enough so that there is no direct 
tunneling but close enough so that the Coulomb interac- 
tion between the surfaces is not negligible. In this case, 
the Majorana signal can be used to infer the existence of 
interaction-mediated tunneling between the surfaces. 

In the sandwich structure, when there is direct tun- 
neling between the two surfaces, bulk TI behavior is ex- 
pected to permeate across the spacer. There is a Ma- 
jorana signal in this situation. Nevertheless, it is possi- 
ble to differentiate this single-particle tunneling from the 
interaction-mediated tunneling of TEC by studying the 
temperature dependence of the Majorana signal. The im- 
portant point is that the bulk TI gap, which is ^ 300 meV 
in Bi2Se3, is much larger than the typical TEC gap, which 
is expected to be '-^ 0.1 meV. Therefore, a Majorana sig- 
nal caused by the single-particle tunneling would be un- 



usually robust as the temperature is increased beyond 
the TEC gap. It should, however, be pointed out that it 
is unlikely that the thickness necessary to allow separate 
superconducting contacts inside the spacer will be in this 
range. 
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Gilbert, and T. L. Hughes. I thank H. Fertig and S. 
Vishveshwara for useful comments and G. Ortiz for the 
same and also for bringing Ref. [23]to my attention. This 
research was supported by the College of Arts and Sci- 
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Note added. — In the final stages of this project, a 
preprint by Meng, Vishveshwara and Hughes |25j ap- 
peared where similar mass terms have been considered. 
However, they focus on the boundary between a magnetic 
domain and a tunneling region. 
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